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Spin dependent single electron tunneling in ferromagnetic double junc- 
tions is analysed theoretically in the limit of sequential tunneling. The in- 
fluence of discrete energy spectrum of the central electrode (island) on the 
spin accumulation, spin fluctuations and tunnel magnetoresistance is anal- 
ysed numerically in the case of a nonmagnetic island. It is shown that spin 
fluctuations are significant in magnetic as well as in nonmagnetic junctions. 



Single electron tunneling in mesoscopic double-junctions has been recently 
extensively studied both experimentally and theoretically. 1-3 For a small ca- 
pacitance C of the central electrode (island) the charging energy E c = e 2 /2C 
can be larger than the thermal energy ksT. Discrete charging of the island 
leads then to Coulomb blockade of the electric current at voltages below a cer- 
tain threshold voltage, and to a characteristic 'Coulomb staircase' at higher 
voltages. When the number of electrons on the island is not large, the energy 
quantization may become important as well, and can lead to additional steps 
in the I — V curves (or additional peaks in the dl/dV characteristics). 4 ' 5 

The interplay of ferromagnetism and discrete charging was studied only 
very recently. 6 " 11 Theoretical results show that discrete charging leads to os- 
cillations in tunnel magnetoresistance (TMR) (that is in the change of the 
total junction resistance when relative orientation of the magnetization of 
both electrodes and of the island is varied). 9 ' 11 In Ref.[9] the intrinsic spin 
relaxation time on the island was assumed to be sufficiently short (shorter 
than the time between two successive tunneling events) to neglect spin accu- 
mulation. This model has then been extended to describe the regime where 
intrinsic spin relaxation time on the island is larger than the time between 
successive tunneling events. In this regime spin accumulation has to be taken 
into account, which leads to enhanced TMR and several new phenomena, like 
inverse TMR or negative differential resistance. Spin accumulation can also 
generate TMR when the island is nonmagnetic. 

In the model of Refs.[9,ll] the discrete nature of the density of states on 
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the island and also the fluctuations of the spin accumulation were neglected. 
These restrictions are relaxed in the present paper, where spin fluctuations as 
well as discrete structure of the density of single-electron states are taken into 
account. On the other hand, as in the models refered to above, we restrict 
our considerations to the limit of sequential tunneling (where the orthodox 
tunneling theory is applicable), that is to the limit when the resistance of 
each single junction is larger than the quantum resistance R q , R q = h/e 2 . In 
that limit higher order processes (cotunneling) can be neglected. However, 
these processes may play an important role in the blockade regime. 10 

Geometry of the junction and energy structure is shown schematically in 
Fig.l. In a general case both electrodes and the island can be ferromagnetic. 
We consider the case when the resistances of both left and right junctions 
are much larger than the quantum resistance R q . Moreover, we assume that 
ksT » T, where T denotes width of the discrete energy levels on the island. 
To calculate the electric current / which flows through the junction when a 
voltage V is applied between the source (right) and drain (left) electrodes, 
one can then use the orthodox tunneling theory. 4 ' 5 ' 12-14 

Assume that there are iV* excess electrons on the island. Generally, one 
can write N* = N? + N*, where N* (a =f or a =|) is the number of excess 
electrons of a given spin orientation. The excess spin polarization of the 
island is then equal to N* — N*. If N a is the number of electrons with spin a 
on the island when a voltage V is applied, and N® is this number at V = 0, 
then one can write N* = N a — N®. 
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Denoting by E ia the discrete energy levels of the island, one can express 
the electric current / in the stationary state as 

/ = e EEEI^I 2p (M) 

o i {n} 

x {^ cti1 [1 - f(E ia + Efr - E F )\ - 5 ntafl f{E la + E l +. -E F )} . (1) 

In the above equation P({n}) is the probability of a particular configu- 
ration of the occupation numbers of discrete energy levels of the island, 

{n} = {n 1T ,n 2T , ....,n u ,n 2i , }, with n ia = 1 (n ia = 0) for occupied 

(empty) states and the sum over {n} denotes the summation over all possi- 
ble occupation configurations. This probability can be determined from an 
appropriate master equation. 4,5 Apart from this, T\ a in Eq.(l) is the matrix 
element corresponding to the electron tunneling to (from) the left electrode 
from (to) the level E ia of the island, while E 1 ^* is defined as E 1 ^, = eV^* ±E C . 
Here, V^* = (C r /C)V+N*e/C, where C r is the capacitance of the right junc- 
tion, C is the total capacitance of the island, and — e is the electron charge 
(e > 0). Finally, we assumed in Eq.(l) the equilibrium Fermi-Dirac distribu- 
tion of electrons in the electrodes, with the Fermi level Ep at vanishing volt- 
age. In the following, we will also introduce the quantity E r ^i = eV^* T E c , 
where Vfc, = (Ci/C)V — N*e/C and C\ is the capacitance of the left junction 
(C = C\ + C r ), and the matrix element T[ a for tunneling through the right 
junction. 

It is convenient to introduce the probability P(N^,N^) of finding on the 
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island TVj electrons with spin a —\ and N± electrons with spin a =[, 

P(N hNl ) = E^W^fEp-pt <WE P «* > ( 2 ) 
{«} 

and the distribution function 

F(E ta \N hNl ) = pT^rrE^W)^ <We,«,t <^E^r ( 3 ) 

V I ) 1/ { n } 

Equation (3) describes the probability that the level E^ is occupied under 
the condition that the island contains jVj electrons with spin a —\ and iVj_ 
electrons with a —\. The current / is then given by 

/ = - e EEEE^^){[ 1 -WT.^)]W(^ + 4-^) 

- P(P i(J |iV T) JV A )|!Z£| 2 [1 " f(E tu + E 1 ^ - E F )}} . (4) 

In the stationary state one finds the following master equation for the 
probability P(N^, JVjJ: 

^P(iV T , iVj = = -P(N h 7VjA(iV T , iVj 

+P(iV T + 1, A^P^ + 1, Ni) + P(N h N t + 1)5^, ^ + 1) 
+ P(iV T - 1, JVJC^ - 1, N t ) + P(iV T , N t - l)Q(iV T , ^ - 1) + ft s , (5) 
where 

MN h N t ) = E Et 1 - n^l^ T , ^)]{|7tlV(^ + - ^) 

<r i 

<y i 
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+ \T[ a \ 2 [l-f(E la -E r N --E F )}} 



(6) 



Ber(Nf, iVj_) 



E *W#T» ^){l^| 2 [l - + E l m - E F )\ 



+ \T[ a \ 2 [l-f(E ia -E r N --E F )}} 



(7) 



C CT (7V T ,^) 



£[1 " ^)]{|^| 2 /(^ + - ^) 



+ \T[ a \ 2 f(E ia - E& - E F )} . 



(8) 



In Eq.(5) 1Z S stands for terms responsible for magnetic relaxation on the 
island. Note that spin- conserving relaxation processes are included in Eq.(5) 
only through the distribution function F{E ia \N^Ni). 

We assume in the following that the energy relaxation time due to spin- 
conserving relaxation processes is significantly shorter than the time between 
successive tunneling events, and also shorter than the spin relaxation time. 
This allows us to use an equilibrium form F^E^N^N^) of the distribu- 
tion function F{E i(T \N^Ni), which depends on the spin orientation. In a 
general case the equilibrium distribution F eCL (E i(T \N^Ni) can be found from 
the Gibbs distribution. In the limit k B T 3> AE the distribution function 
Feq^io- 1 iVfiVj) is equal to by the Fermi-Dirac distribution 



where the chemical potential /i a (N a ) is determined by the equation 



As noticed by Beenakker [5], the equilibrium distribution for k B T ps AE 
differs significantly from the Fermi distribution. 



F{E la \N h N i )=f(E la - l i a {N a )) , 



(9) 



Y t f{E iff - l i a {N a )) = N a . 



(10) 
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Consider now some numerical results. For simplicity, we will restrict the 
following considerations to the situation when the island is nonmagnetic and 
the magnetizations of the electrodes can be either parallel or antiparallel. 
To emphasize the role of spin accumulation we assume that the intrinsic 
spin relaxation time on the island is long enough to neglect all intrinsic 
spin-flip processes. The magnetic relaxation on the island takes then place 
only through tunneling processes. For simplicity we assume that the energy 
levels on the island are equidistinct, with the inter-level spacing AE. Apart 
from this, we assume that ksT < AE (but still ksT » Y) and AE < 
E c . In Fig. 2a we show the I-V characteristics for parallel and antiparallel 
configurations. In both cases the electric current is blocked below a threshold 
voltage equal approximately to 13 mV. Above the threshold voltage typical 
'Coulomb staircase' appears with additional small steps due to discrete levels. 
The existence of those steps leads to additional peaks in the derivative dl/dV, 
as shown in Fig. 3b for the antiparallel configuration (the corresponding curve 
for the parallel configuration was not shown there for clarity). 

The difference between I-V curves for the parallel and antiparallel con- 
figurations is due to a different spin accumulation in both geometries (note, 
that both I-V curves become identical when no spin accumulation occurs in 
the nonmagnetic island 9 ' 11 ). To show correlations between the I-V curves 
and spin accumulation, we present in Fig. 2c the average value of the differ- 
ence between the numbers of spin-up and spin-down excess electrons on the 
island, (iV? — N*), for both antiparallel and parallel configurations (up to a 
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constant factor, Fig. 2c shows an average magnetic moment induced on the 
island due to spin accumulation). For the symmetrical structure assumed 
here (R^/R^ = R r -\/R r i for the parallel configuration, where Ri( r ) a are the 
junction resistances for a given spin channel) there is no significant spin ac- 
cumulation in the parallel configuration. Only small nonzero values occur 
around the Coulomb steps. For asymmetrical junctions spin accumulation 
also occurs in the parallel configuration. Discrete structure of the density of 
states on the island appears in the voltage dependence of (TV* — TV*) as the 
fine steps clearly visible in Fig. 2c. 

The number N* — N*, of spins accumulated on the island fluctuates 
around its average value {N* — TV*), as shown in Fig. 2d, where the stan- 
dard deviation [((iV* — N*) 2 ) — (N* — TV] 1 ) 2 ] 1 / 2 is plotted against the voltage 
V. It is worth to note that although there is almost no spin accumulation in 
the parallel configuration, the corresponding fluctuations are relatively large 
and comparable with the fluctuations in the antiparallel configuration. The 
fine structure in the voltage dependence of the standard deviation originates 
from the discrete electronic structure of the island. 

As we have already pointed above, spin accumulation on the island gives 
rise to a difference between the I-V curves in the parallel and antiparallel 
configurations. This difference, in turn, leads to the tunnel magnetoresis- 
tance (TMR) defined quantitatively as (R ap — R p )/R p , where R ap and R p 
are the total junction resistances respectively in the antiparallel and parallel 
configurations. The ratio (R ap — R p )/R p is shown in Fig.2e, where the broad 
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peaks correspond to the Coulomb steps while the fine structure originates 
from the discrete structure of the density of states. 

In the limit of a nonmagnetic junction (nonmagnetic island and nonmag- 
netic electrodes) there is no spin accumulation on the island and no TMR. 
However, spin fluctuations still occur as shown in Fig. 3 for the limit of no 
intrinsic spin relaxation on the island. As in Fig. 2, the discrete energy levels 
on the island lead to fine steps in the voltage dependence of the standard 
deviation. 

In summary, we developed in this letter a formalism for calculating elec- 
tric current and tunnel magnetoresistance in ferromagnetic double junctions, 
which takes into account spin accumulation on the island and discrete struc- 
ture of the density of states. The discrete energy levels lead to fine steps 
in the I — V curves and to additional peaks in the first derivative dl/dV. 
Moreover, we showed that the discrete levels give rise to fine structure in the 
spin accumulation and spin fluctuations on the island, as well as in TMR. 
To single out the charging effects we neglected in this letter tunneling pro- 
cesses with simultaneous spin flip, which lower TMR and also lead to a smoth 
decrease of TMR with increasing voltage. 15 ' 16 

The work was supported by the Polish Research Committee through the 
research project 2 P03B 075 14. 
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1 Schematic of the double-junction considered in this paper. Magnetization 
of the island and of the right electrode can be either parallel or antiparallel 
to the magnetization of the left electrode, as indicated by the solid and open 
arrows. Discrete levels of the island are marked schematically by the solid 
(dashed) lines for spin-up (spin-down) electron states. 
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Fig. 2 Voltage dependence of the tunnel current / (a), first derivative dl/dV (b), 
spin accumulation (N* - N*) (c), standard deviation [((N* - N^) 2 ) - (N* - 
iV*) 2 ] 1/2 ( d ); and tunnel magnetoresistance TMR (e), calculated at T = 2.3K. 
The solid and dotted curves in (a), (c) and (d) corespond to the antiparallel 
and parallel configurations, respectively. The other parameters assumed in 
the numerical calculations are: AE = 3meV, C\ = 6.6aF, C<i = 1.32aF 
(which gives E c = lO.lmeV), = 200Mfi, R n = lOOMfi, R rl = 4Mfi and 
i? rT = 2Mf2 for the parallel configuration (R rl = 2Mfi, i? rT = 4Mf2 for the 
antiparallel configuration). ^ 
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Fig. 3 Fluctuations of the spin accumulated on the island in the limit of a non- 
magnetic double junction. Standard deviation [((N^ — N*) 2 ) - (N^ - N*) 2 ] 1/2 
is shown as a function of the voltage V for the parameters: T = 2.3K, AE = 
3meV, C x = 6.6aF, C 2 = 1.32aF (E c = lO.lmeV), R n = R n = lOOMfi, and 
R r1 = R, = 4Mfi. 
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